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Abstract. We obtain Lipschitz regularity results for a fairly general 
class of nonlinear first-order PDEs. These equations arise from the inner 
variation of certain energy integrals. Even in the simplest model case of 
the Dirichlet energy the inner- stationary solutions need not be differ- 
entiable everywhere; the Lipschitz continuity is the best possible. But 
the proofs, even in the Dirichlet case, turn out to relay on topological 
arguments. The appeal to the inner-stationary solutions in this context 
is motivated by the classical problems of existence and regularity of the 
energy-minimal deformations in the theory of harmonic mappings and 
certain mathematical models of nonlinear elasticity; specifically, neo- 
Hookian type problems. 



1. Introduction 

We establish Lipschitz regularity of solutions of nonlinear first-order PDEs 
that arise from inner variation of numerous energy integrals. This includes 
the model case of the Dirichlet energy for mappings /i: — )• $7* between 
two designated domains in C. Roughly speaking, the inner variation of h 
amounts to composing h with a diffeomorphism of Vt onto itself. This type 
of variation is often used when the standard first variation is not allowed. 
For instance, when dealing with mappings with nonnegative Jacobian the 
inner variation is necessary to preserve the sign of the Jacobian. Let us 
begin with the Dirichlet integral. 



(LI) S,,[h]= IJj^^\' = ^ jj^ 



|^2p + l/isP) dxi dx2 , z = xi + ix2 



Hereafter hz = ^ and = ^ are complex partial derivatives of h . The 
first variation of S" results in the Euler-Lagrange equation, 

(1.2) Ah = Ahzz = 
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In contrast, the inner variation leads to a nonlinear equation 

d_ 
dz 



Q 

(1.3) — (hzfi^ = 0, equivalently, hz hz = </>, {4> is analytic) 



for mapping in the Sobolev space '^^^{^) ■ This equation will henceforth be 
referred to as the Hopf-Laplace equation. There are important nonharmonic 
solutions of (1.3). Such solutions arise typically as weak limits of the energy- 
minimizing sequences of diffeomorphisms h: O Passing to the limit 
we often loose harmonicity; at the points where the limit map fails to be 
injective [16], and only at those points [17]. 

The unavailability of the Euler-Lagrange equation is a major source of 
difficulties in the theory of nonlinear elasticity [7, 8, 29]. This drives one 
to investigate the regularity of the energy-minimal mappings on the basis 
of the inner-variational equation alone; also known as energy-momentum or 
equilibrium equations, etc [11, 27, 30]. Several results in this direction were 
obtained in [9, 12, 25, 32]. Nevertheless this theory is still in its infancy. 

The most general setting we are dealing with can be described as follows. 
Let H = T-L{z, ^) be a continuous function in x C^^ , where Cr = : -R < 
1^1 ^ oo }, ^ < oo. We impose two structural conditions on %. The 
first one is the Lipschitz condition with respect to the reciprocal of the ^ 
variable, 

1 _ 1 

6 ~6 



(1.4) \H{z,^{)-n{z,^2)\^L- 



for 2; G and ^1, ^2 G Cr. 



with a constant ^ L < 00. Regarding regularity with respect to 2; € il, we 
shall require that the function z^%{z,^) be Holder continuous. Precisely, 
the second structural condition reads as: 

(1.5) su^\H{z,^)\ + sup 1^(^1^0-^(^2,61 ^j^^ for|e|>i? 

zm zx+zi \Z\ — Z2\°' 

Definition 1.1. A function h G "^iJ^iS^ is said to be a solution of the 
equation 

(1.6) h-z = n{z,hz) 

if (1.6) holds for almost every point z G O , whenever |/t2(2;)| > R- 

Note we impose no condition at the points where |/i2(2;)| ^ R; this yields 
boundedness of the gradient of /t, [/i^l ^ [/i^l ^ ^ . 

A simplified version of our main result reads as follows 

Theorem 1.2. Let the equation (1.6) comply with the conditions (1.4) 
and (1.5). Then every solution h G W^''^{Q,) with nonnegative Jacobian 
is locally Lipschitz continuous but not necessarily "i^^ -smooth. 

A special case of Theorem 1.2 deserves a separate mention because it 
covers the variant of (1.3) with not necessarily analytic right-hand side <p. 
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Theorem 1.3. Let h G #^^'^(0) be a mapping with nonnegative Jacobian. 
Suppose that the Hopf product hz hz is bounded and Holder continuous. Then 
h is locally Lipschitz but not necessarily -smooth. 

More specific statements, including gradient estimates near , are pre- 
sented in Theorem 3.1. Examples of variational problems leading to equa- 
tions in Theorem 1.2 are provided in Section 2. 

Our proofs draw upon the theory of Beltrami equations combined with 
methods of topology, a technique originated in the theory of general nonlin- 
ear first order elliptic systems [3, 5, 10, 15]. Initially, for the linear elliptic 
system 

hz = nhz + v hz, ImI -\-\v\^k <\ 
the solutions are quasiregular, a concept firmly rooted in the geometry of 
analytic functions. While it might sound trivial, the observation that the 
difference of two solutions is again quasiregular is deep and useful, for it pro- 
vides us with powerful topological tools to obtain existence, uniqueness and 
regularity of solutions. And this was exactly a blueprint for the definition 
of ellipticity of fully nonlinear Beltrami type equations, 

(1.7) h-z = n{z,h,), \n{z,^i)-n{z,^2)\^mi-^2\, o^A;<i 

Because of nonlinearity the difference of two solutions need not solve the 
same equation but it does solve another elliptic equation and, as such, sat- 
isfies the distortion inequality 

(1-8) \h-9-z\ ^ k\fz-gz\ 

meaning that f — g \s quasiregular, see Definition 3.2. 

A chief distinction from the elliptic cases discussed above is that the 
solutions of (1.6) need not be quasiregular. However, we will construct a 
continuous family {-F^}agc of "good" solutions of (1.6) such that — 
h are quasiregular. After that we appeal to the topological properties of 
quasiregular mappings. 

We believe that these ideas will have applications to even more general 
PDEs than those in our paper. An interested reader is referred to recent 
papers by D. Faraco, B. Kirchheim and L. Szckclyhidi [13, 21] which also 
combine the theory of quasiregular mappings with topological arguments. 

2. Inner-variational equations 
Let us consider the energy integral for mappings h :Q, ^ C 

(2.1) = // ^{z,h,hz,hz) dx\dx2 , z = x\-\-ix2 

J Jn 

where E = E{z,w,^,C) is a given stored-energy function. Prom the point 
of view of Geometric Function Theory the mappings h must take J7 onto a 
designated domain Q* . In the elasticity theory these domains are referred 
to as the reference configuration and deformed configuration, respectively. In 
the neo-Hookean model of hyperelasticity the stored-energy function blows 
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up when the Jacobian of h approaches zero. Thus one is looking for mapping 
h with positive Jacobian determinant which minimizes the energy. It is not 
always the case that the minimizers satisfy the Euler-Lagrange equation. 
However, it is legitimate to perform the inner variation of the energy integral. 

Given any test function rj G '^^{^l) and a complex parameter t, small 
enough so that the map z z + trj{z) represents a diffcomorphism of O 
onto itself, consider the inner variation h^{z) = h{z + trj) and its energy 



<r[/t*] = B{z,h\hl,hl)dxidx2 
J Jn 



First we make a substitution w = z + trj{z) and then differentiate to 
obtain an integral form of the equilibrium equation ■§i\f^Q'S'[h^] = . We 
eliminate rj through integration by parts to arrive at what is called the 
inner-variational equation 



d 



d 



(2.2) — h,B^ + /i^E^- + — h,B^ + /i^E^- - E 



dz 



+ E, =0 



Here the subscripts under E stand for complex partial derivatives of E = 
'Ei{z,w,^X) ■ The partial derivatives ^ and ^ are understood in the sense 
of distributions. Note that does not enter this equation; meaning that 
E is not required to be differentiable with respect to the u;-variable. We 
say that h is inner stationary for S if it satisfies the equation (2.2), 
regardless of whether h is energy-extremal or not. 

2.1. Dirichlet integral. The most basic example is the Dirichlet integrand 
E = l^p + l^p and the associated Hopf-Laplace equation (1.3), 

d_ 

dz 



{h,h,)=0, for heW,l'^{n) 



2.2. Poincare disk. More general equations arise in the theory of harmonic 
mappings between Riemann surfaces [19, 20, 23, 28]. In particular, let the 
target be the Poincare disk. This is the unit disk J]) = {w e C: \w\ < 1} 
equipped with the hyperbolic metric ds = j^q^- The associated Dirichlet 
integral 



J Jo, 



dxi dx2 



In 

is certainly infinite for homeomorphisms h : O D in the Sobolev space 
■^oc^(^'^) • Nonetheless, it is interesting to examine the inner-variational 
equation and all its solutions, not necessarily homeomorphisms. 

For slightly greater generality we consider the weighted Dirichlet integral 
m= II {\h,\^ + \h,f)p{z,h)dxidx2 
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and its inner-variational equation 

(2.3) —[piz,h)h,Y,] = p,iz,h){\h,\^ + \h,\''), ioi heW,l'^in,]D)) 



Theorem 2.1. Suppose p = p{z,w) ^ 1 is Lipschitz continuous in the 
z-variable and Holder continuous in the w-variable. If h & ^\o'c i^^^) ^-^ 
a solution of (2.3) with nonnegative Jacobian, then h is locally Lipschitz 
continuous. 

Proof. Generally speaking we are dealing with a nonhomogeneous Cauchy- 
Riemann equation 

QU 

(2.4) -Q^ = u, where U = p{z,h) h^h^ , u = {\h^\'^ + Ih^l"^) p^{z,h) 

The equation is elliptic, so we gain some regularity properties of U and u. 
At the beginning we only know that U,u E ^^^{Q,). We shall recurrently 
improve integrability properties of these terms. First observe that U, having 
^-derivative in ^^^^(Q), lies in ^^^^{^) for every exponent 1 < s < 2. 
Then, in view of pointwise inequality ^ p{z,h) \ hz\\hz\ = \U\, we see 

that G -S?ioc(^) • This implies that also € ^il^i^) ■ In this 

way we gain higher integrability of the right hand side of (2.4); namely, 
u = {\hz\'^ + l^sP) Pzi^jh) G ^[^^(ri) , because pz{z,h) is bounded. Now 

2s 

equation (2.4) places U in the space ^i^^' {^) ■ This, in view of P ^ l^^li 

23 23 

yields \hg\'^ G ^Sfj^-" (Q) . Hence G ^^^-^ (n) as well. Thus we gained 

more integrability of u; namely, u G ^^^{^) , with p = > 2 . We again 
turn back to equation (2.4). This time the equation yields Holder continuity 
of U ; that is, U G '^^^i^) with a = 1 — | > . Let us write the equation 
as 

p{z, h) 

We observe that h is also locally Holder continuous, because hz G ^^J^{VL) 
with exponent 2p > 2. The conclusion is that the Hopf product hz hz is a 
Holder continuous function. By Theorem 1.3, h is locally Lipschitz. □ 

2.3. An application to nonlinear elasticity. We now turn to some fairly 
general energy integrals of interest in nonlinear elasticity. Given two bounded 

domains Vl and J7* in C, we consider mappings h: O — )■ of Sobolev class 

W^''^{^) whose Jacobian determinant Jh == detDh = \hz\^ — \hz\^ is non- 
negative. In nonlinear elasticity of isotropic materials one considers the 
energy of h of the form 

s[h] = [[ w{z,h,\hz\Mh-z\^) 

J Jn 

Specifically, neo-Hookean models of elasticity [6] deal with the integrands 
W which blow up as the Jacobian determinant approaches zero. For the 
sake of simplicity we forgo the dependence of W on the z and h variables. 
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The interested reader may generalize our considerations by including z and 
h variables to the integrand, like in §2.2. For the energy integrand 

E(m) = t^(l^,lM^,l') 

the inner-variational equation (2.2) simplifies as follows 

(2.5) ^ [(^« + ^b)h,h;] + ^ [\h,\^Wa + \h-,\^Wb -W]=Q 

where W and its partial derivatives Wa and are evaluated &t a = \hz\'^ 
and h = \hz\^. To emphasize a possible neo-Hookean character of the inte- 
grand we express it as 

^ ' ' (a-6)P-i 

so 

where F = F{a,b) is defined and continuous in O = {(a, 6): a ^ 6 ^ 0} 
-the closure of the first octant O = {(a, 6) : a > 6 > 0} . Let us remove the 
corner of O to introduce Oo = O \ {(0, 0)} . As usual, we write 

\VF\ = \Fa\ + \Fb\ and \V^F\ = \Faa\ + \Fab\ + \Fi,i,\ 

We make the following standing assumptions on F : O ^ [0, oo) : 

(2.7) Fg '^(O)n^2(O\{(0,0)}) 

(2.8) F is homogeneous of degree p ; that is, 

F{ta, tb) = tPF{a, b) for t ^ , a^b^O. 

Furthermore, for a > 6 > , we assume that 

(2.9) ^a + br<F<ia + br 

(2.10) (a + bf-^ <Fa + F,^ \\7F\ <{a + bf'^ 

(2.11) IV^FI < (a + 6)f-2 

Here the notation < refers to an inequality with the implied constant (pos- 
itive) which stays independent of (a, b) G O. 

Theorem 2.2. Let h G ^^'^ {^) be an inner-stationary mapping for the 

energy integral (2.6) with S[h\ < oo, where F satisfies the conditions (2.7)- 

(2.11) . Then h is locally Lipschitz continuous. Furthermore E(Dh) is locally 
bounded. 

Proof. First observe that 

(2.12) / / \Dhf < ^[h] < oo 

J Jn 
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hence h G Regarding inner variation, we apply formula (2.5) 

to W{a,b) = (^a-hyp-^ ^^^^ oVFa + bW^ — W = 0, because of p- 

homogeneity of F . The equation (2.5) takes the form 



(2.13) 

where we note that 



d_ 



Fa+F^ 



hzhz 







(a - b)P-^ 



< 



(a - 6)P-i 
(a + b)P-^Vab ^ (a + b)P 



< 



F{a,b) 



{a - b)P-^ ^ (a - b)P-^ ~ (o - b)P-^ 
Thus, by Weyl's lemma, (2.13) yields that 

(2.14) — ° h)P~^ ^zhz = is an analytic function, actually in =Sf ^(f]) 
We also note at this point that in view of condition (2.10) 



(2.15) 



Fa + Fb 



(a + b)P-^ 



Fa + Fb 



(a - 6)P-i 



hzhi 



Hence hg G ^i^(l^). 

We are going to solve (2.14) for hz in terms of (p and hz, at least when 
\hz\ is sufficiently large. Choose and fix an arbitrary subdomain Q' Q and 
let 

= Il0lb-(n') < oo 

Since Fa+Fj, is homogeneous of degree p—l, equation (2.14) can be written 
as 

F„(l,fc2)+F6(l,fc2)7^ 



(2.16) 



(1 - fe2)P-l /l^ 

where ^ = ^ V2 , provided 

(2.17) \hz\ ^ 2N. 

Taking the absolute values of both sides in the equation (2.16), we obtain 

.o1c^ F,(1,P)+F,(1,P) , _ 101 ^ 

(l-P)f-i ''-Jf^-' 

The left hand side represents a "^^-smooth function $ = $(A;), < A; ^ 1/2 . 
We see that 

$'(0) = F,(1,0) + F5(1,0)>1 

by the condition (2.10). Therefore (2.18) admits unique solution for k close 

to 0, say 

k = sT{s), 
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provided ^ k ^ ko ^ V2, where T = r(s) is a "^^-function in a small 
interval ^ s ^ So- Now the equation (2.16) reads as, 

(2.19) V, = ^r(^), whenever ^ £ 

Thus we arrive at the equation of the form (1.6), where 



All the conditions on H in Theorem 1.2 are satisfied, so 

\Vh\ G ifi-(^7'). 
It remains to estimate the integrand. We have the identity 
B{Df) _ F(l,fc2) 



\<P\ + \h,\ (1 - P)P-1 + [Fail, k^) + F,{1, fc2)] k 

regardless of whenever \hz\ ^ R holds or not, where k = |^ € [0,1]. The 
right hand side is bounded. Therefore, 

B{Dh) < \cp\ + \h,\ G □ 

2.3.1. An example. The class of energies in (2.6) covers the following par- 
ticular integral, 

. ,., rr \Dhp rr (ift.i' + i^i')" _ , 

subject to homeomorphisms h: in the Sobolev space W^''^{Q,). 

This case gains additional interest in Geometric Function Theory because 
the transition to the energy of the inverse mapping / = : O* 
results in the ^^'-norm of the distortion function 

We see that conformal mappings, for which Kf = 1 , are the absolute mini- 
mizcrs. In general, .jSf^-integrability of the distortion function only guaran- 
tees that / G Indeed, 



V 



= \\Kfr -Mi^ <oo 

We do not pursue this matter further; see [2, 4, 18, 24] for more on the 
minimization of . 
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2.3.2. Nonisotropic energies. The methods presented in this paper are also 
pertinent to some nonisotropic energies (2.1). By way of illustration, here 
is an example of such energy integral 

where F = F{a, b, c) is a function of three real variables a ^ 6 ^ and 
c G M . We assume, as in §2.3, that F is homogeneous of degree p ; that is, 

F{ta, tb, tc) = tPF{a, b, c) 

An elementary but tedious computation reveals that the inner-variational 
equation (2.2) takes the form : 

{Fa + Fb)hJ~, + {hl + h^^)F, 



- — 1 = 6 -an analytic function in Q, . 

{\h.\^-\h,\r-' 

We leave it to the interested reader to complete this investigation by im- 
posing precise, fairly minimal, conditions on F = F{a, b, c) in order to 
implement Theorem 1.2. The conclusion is that h is locally Lipschitz and 
the integrand E(D/i) is locally bounded. An explicit example is: 

F{a,b,c) = {a + b)P + €c'^{a + b)P-'^, p^2 

with e > sufficiently small. 

3. Elaborate statement and examples 

The nonnegativity of Jacobian for a solution of (1.6), under the structural 
conditions (1.4)-(1.5), implies 

(3.1) II/12II00 < M + \/Z + i? 

For the sake of greater generality, and clarity of the arguments as well, we 
reformulate our main result (Theorem 1.2) replacing the Jacobian condition 
with the boundcdncss of hg. Furthermore, we give a quantitative result with 
the sharp asymptotic bound on the gradient of h near the boundary. Define 



and 



OSCj, [h] = sup I h{a) - h{b) 



\Vhiz)\^\h,iz)\ + \h,{z)\ 



Theorem 3.1. To every equation (1.6) there corresponds a structural con- 
stant Ao = Ao(H) such that the following holds. Suppose that a function 
h e W^''^{n)n^^{n) with G ^°°in) satisfies (1.6). Then h is locally 
Lipschitz. Moreover, for almost every z & Q we have 

(3.2) \^Kz)\ ^^osc^[h] + 4||/i5||^oo(n) + 6A0 



10 T. IWANIEC, L. V. KOVALEV, AND J. ONNINEN 

where r = min{dist(2;, dQ), 1}. In particular, 

limsup \Vh{z)\ ■ dist(z, 90) ^ 3 OSC^[h] 
2-> an 

If, in addition, h is continuous up to the boundary then 

limsup \Vh{z)\ ■dist{z,dn) =0 
z-^dQ 

Prerequisites. Before proceeding to the proof of Theorem 3.1 we need some 
definitions. 

Definition 3.2. A mapping g G W^^'^{VL) is quasiregular if there exists a 
constant k < 1 such that \gz\ ^ k\gz\ a.e. in O. Such a mapping is cahed 
quasiconformal if it is also injective. 

The Jacobian determinant Jg = \gz\'^ — [gzl"^ of a quasiregular mapping 
is clearly nonnegative. It is well-known that a nonconstant quasiregular 
mapping is both open and discrete, sec [3, Corollary 5.5.2] or [22, Chapter 
VI]. In the sequel we shall appeal to the following topological fact which 
relates the cardinality of preimages with the topological degree [14, Theorem 
4] and [26, Proposition 1.4.10]. 

Lemma 3.3. Letil G C be a bounded domain and G: CI C an orientation 
preserving continuous mapping that is open and discrete in Jl. Then for 
every v ^ G{dCl) we have 

Card{z G n -. G{z) = v} ^ degJv-G] 

Examples. Here we provide examples that demonstrate the sharpness of 
both assumptions and conclusions of our theorems. 

The first example shows that even the most basic of our equations, 

hzhz = -1 

admits /^-"^'^ -solutions with ^ that are not "^^-smooth. 
Example 3.4. Let n = {z: \z-l\ < 1/2}. For 2; G O \ [1,3/2) we define 

h{z) = - z'/') + (1 - z'/'f' + (1 - f/'f' 

using the principal branch of the power function. Then 

Therefore, hzhz = —1- Note h has a Lipschitz extension to Cl, in fact h = 
on the interval [1,3/2)- One can also check that ^ a.e. However, 
h ^ ^^{O,), because z^/^{hz + hz) = 2(1 — z^l^'fl'^ fails to be continuous. 
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A much simpler example can be given if one does not insist on the special 
Hopf-Laplace structure hzhz , arising from the Dirichlet energy integral. Let 



Hz) 




Imz ^ 
Imz < 



This is a bi-Lipschitz mapping of C which is not "^^-smooth but satisfies the 
equation 

hz = j-- 2, 

as well as any equation of the form hz = T-L{hz) with H(2) = 1 and H(3) = 0. 
Theorem 3.1 also encompasses other Hopf type products such as 

\hz\h, = 4^, i^e^'^m, 

but fails for some Hopf type products, like in the pseudo Hopf-Laplace equa- 
tion below. 

Example 3.5. The equation 

(3.4) h,\hz\ = 1 

admits a solution that is quasiconformal but not locally Lipschitz. 
Proof. We look for h in the form 

(3.5) h{z) = 2ztP{-2\og\z\), \z\ ^ 1 

where ijj: [0, oo) [l,oo) is a strictly increasing function with ■0(0) = 1. 
Since 

2z 

(3.6) /i, = 2V(-21og|z|)-2V''(-21og|z|) = -_V;'(-21ogH) 
the function must satisfy the differential equation 

(3.7) ^^^t)-^P'{t))^P'{t) = ^, < t < oo. 
Rewriting (3.7) as an equation for the inverse function 

(3.8) |; = ^(^ + ^^^)' *^i) = °' 

we arrive at 



(3.9) t = ip^ -1 + V'VV'^ - 1 - log(V' + VV'^ - 1). 

The equation (3.9) determines a differentiable strictly increasing function 
[0,oo) [l,oo), since the right hand side of (3.8) is positive. Note that 
ip(t) — >■ oo as t — >■ oo, hence h is not Lipschitz in any neighborhood of the 
origin. In view of equation (3.7) this implies tp'^t) — ^ as t — ?■ oo. It now 
follows from (3.6) that h is quasiconformal in a neighborhood of the origin; 
in particular Jh > almost everywhere. □ 
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The major difference between equations (3.3) and (3.4) is that the latter 
is not solvable for hz in terms of hz ■ 

Our final example shows that the Holder continuity of z i— t- Ti^z,^) in 
Theorem 3.1 cannot be relaxed to continuity, even for the standard Hopf 
product. 

Example 3.6. Let h{z) = zloglog|z|~^ for \z\ < 1/2. This mapping is an 
orientation preserving homeomorphism which belongs to W^'^ for all p < 00. 
We compute 

11 z 1 

hz = log log —2 - log"^ —2 and hz = - log"^ 

z Izr 



Clearly, hzhz is continuous. However, h is not Lipschitz. 

Even for the most basic equation h^h^ = 1, allowing the Jacobian of h to 
change sign destroys any hope for improved regularity [12]. 

4. Model case: the Hopf-Laplace equation 

In order to illustrate our ideas without getting into technicalities, we first 
take on stage the Hopf-Laplace equation 

(4.1) hz hz = (f){z), with analytic in O C C 

Theorem 4.1. Suppose that the Hopf product hzhz = (l>{z) is analytic 
and bounded in a domain Q C C , for some h G n with 

hz G =Sf°°(ri). Then h is locally Lipschitz. Moreover, for almost every 
z E we have 

^'•'^ l^^^^^l ^ ZZ^dn] + 2||/..|k-(a) + 3||<^||^i(^) 

Proof. Finding good solutions to equation (4.1) presents no difficulty. First 
consider = D -the unit disk, and assume that (j) is bounded in B . 
Denote by $ = ^{z) the anti-derivative of (p such that ^(0) = 0. Thus 
$2 = and $2 = ^ . Clearly, $ extends continuously to the closed unit 
disk D = B. The mappings F^{z) = Xz + f^{z) , where f^{z) = X-^^{z) 
with complex parameter X =^ , solve the same equation (4.1). Also note 
that ||/'*'||oo ^ |-^|~^||^!'||oo • A short computation reveals that the difference 
g = g^{z) = F^{z) — h{z) is a >^^'^(B)-solution to a linear Beltrami type 
equation 



gz{z) = u{z) gz{z) , u{z) = — = — , \iy{z)\ - 

whenever |A| ^ 2||/i2||oo. Now consider a continuous family of mappings 
G^{z) = i g^{z) = z + j [f^{z) - h{z)] . We have 

\G>^(z) ^kII'^IIoo , ll^lloo ^1 
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provided |A| ^ 2 y^||^||oo and |A| ^ 13 ||/i||oo • This shows, in particular, 
that is a nonconstant quasiregular mapping, thus orientation-preserving, 
open and discrete. At this point we appeal to a Rouche type lemma. 

Lemma 4.2. Let G = G^{z) be a continuous family of mappings G^ : 
D — )• C parametrized by complex numbers X with g ^\X\ ^oo , such that 

(i) G^{z) = z 

(ii) \G^{z) -z\<l , for zeG and |A| ^ g 

(iii) For every |A| ^ the map G^ : E) — >■ C is orientation preserving 
open and discrete. 

Then, given any € and parameter |A| ^ g, the equation 

(4.3) G^{z) = G^{zo) , /or z G D ; 
admits exactly one solution z = Zq . 

Proof of Lemma 4.2. For \z\ = 1, we see from (ii) that |G^(2;)| > |. Fix 
any point v of modulus |w| < | , so ^ G^(5D) . Therefore, we have well 
defined topological degree of denoted by 

degjt;; G^] for |A| ^ ^> 

This is an integer-valued function, continuous in A , thus constant. For 
A = DO the degree is 1, because G°°{z) = z . Therefore 

deg^[i;; G'^] = 1 , for all parameters |A| ^ ^. 

Now comes another topological fact concerning orientation preserving 
open discrete maps. It asserts that the cardinality of the preimage of v 
does not exceed the degree, see Lemma 3.3. In symbols, 

^ Card{z G D: G^{z) = v} ^ deg^b; G^] = 1 , for all |A| ^ g. 

This applies to the point v = G^{zo) , because \G'^{zo)\ ^ \G^{zo) — Zo\ + 
I -^o I < I + ^ = § , by condition (ii). The lemma is established. □ 

Returning to the proof of Theorem 4.1, we infer that the mappings G^{z) = 
^g^{z) are injective in the disk jD. So are the mappings g^{z) = \z + 
\~^(t{z) — h{z) . This reads as follows: 

(4.4) h{z^) - h{z2) ^ \ .[z^-Z2 + ^ [^{Z^) - $(^2)] } 

for z\ 7^ Z2 in the disk i D . Letting A run over a circle of radius |A| we 
conclude that 

(4.5) \h{zr) - h{z2)\ ^ |A| .\z^-z2 + ^ [$(zi) - $(^2)] 
This is possible only when 

(4.6) \h{zi) - h{z2)\ < |A| .\z,-Z2 + ^ [$(Z1) - $(Z2)] 
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because the right hand side is continuous in A and the inequahty (4.6) holds 
for large values of |A| . 

A conclusion is immediate; the solution h is Lipschitz continuous in the 
disk I D . Moreover, 

(4.7) l|V/i||^oo(iD) ^ 1^1 + ^ U\W^{r>) 

All the conditions we have encountered for the parameter A are satisfied if 
we set 



(4.8) |A| = max 



Therefore, 

(4.9) ||V/i||^oc(iD) ^ 2||/i,-|b-(D) + 13||/t||^^(D) + 3||</)||^i(r,) 

completing the analysis of the case O = B. 

Now let be a general domain. Suppose B(a,r) = {z: \z — a\ ^ r} is a 
closed disk contained in O and /i is a solution to the Hopf-Laplace equation 
in O . We scale down the variables to introduce a function h(z) = ^h{rz+a) 
which satisfies the Hopf-Laplace equation hz hz = (t){rz + a) in the unit 
disk. Inequality (4.9), applied to ^, yields 

13 

(4.10) ||V/l||^oo(iB) ^ 2 ||/t^||_s^c«(B) + y ||/l||^o°(B) + 3 ||(/)||^L(B) 

We are allowed to subtract any constant from h , given that h appears 
in (4.1) only with its derivatives. The estimate (4.2) follows. □ 

Corollary 4.3. Under the hypotheses of Theorem 4.1 we have 

limsup \Vh{z)\ ■ dist(z,(9Q) ^ 13 OSCj,[/i] 

//, in addition, h is continuous up to the boundary then 

(4.11) limsup |V/t(z)| •dist(z,5J^) =0 

The property (4.11) may fail for solutions that are not continuous up to 
the boundary. 

Example 4.4. Here is a bounded solution h{z) = z'^ + sin log z in a half- 
disk 0, = {z: Re2;>0, |z|<l} to the Hopf-Laplace equation 

hzhz = 2 cos log z -analytic and bounded 

This solution exhibits large oscillations arbitrarily close to the point a = G 
dVt. It is for this reason that h fails to have a limit lim^-^o (Rez)|V/i(2;)| . 
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5. Outline of the proof of Theorem 3.1 

The arguments presented for the proof of Theorem 4.1 contain many of 
the key ideas of the proof of our most general result, Theorem 3.1. The 
term 3 ||?!'||^L(n) in (4.2) will be replaced by a number Ao = Ao('H) which 
depends on the equation; that is, the conditions on H. Let us emphasize 
that Ao will not depend on the solution h. 

However, the situation is more intricate because we need to find a coun- 
terpart of the antidcrivativc of (p. The proof of Theorem 4.1 suggests that 
we should look for the family F^{z) = Xz + f'^{z) that complies with the 
equation (1.6); that is, fz = 7i{z,X + fz)- We must show that the latter 
equation admits a continuous family {/^} of "good" solutions. The key 
is that will enjoy uniform Lipschitz bounds, independent of A . Then 
the proof of the Lipschitz regularity of h will be carried out by topological 
arguments in much the same way as in the case of Hopf-Laplace equation. 

We first consider the case Q = D, and treat general domains CI by rescal- 
ing. The sharpness of Lipschitz regularity was already demonstrated in 
Section 3. 



6. Good solutions 

We are looking for a family {/'^: A G C, |A| ^ Ao} of solutions to the 

equation 

(6.1) fz = n{z,x + fz) 

in the closed unit disk D = {z: |2;| ^ 1 }, where H = T-L{z,^) is continuous 
in D X and satisfies the structural conditions (1.4) and (1.5) in J7 = D. 

Proposition 6.1. There is Ao = Xo{a,L,M,R) and a family {/^}|;^|>;^^ 
of solutions to the equation (6.1) in D such that 



(6.2) f\0) = 

(6.3) \fHz,)-f\z2)\^Xo-\zi-Z2\ 

(6.4) |/'H^)-/'^(^)KAo- 

Furthermore, for every z eT) and \X\ ^ Xo , we have 

(6.5) |V/^(z)| ^ Ao 

(6.6) \f^{z)\ , \f^{z)\ ^ l|Ao| ^ ^|A| ^ |A + f^{z) 



Proof. The proof of Proposition 6.1 is proceeded by an extension of equa- 
tion (6.1) to the entire plane C. 
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(6.7) 



H(z, 6 



6.1. Extension to C. We set for all z G C and |^| > R, 

^'H{z,0 if \z\ ^ 1 

{2-\z\) -mi/z,^) ifl^|z|^2 
for \z\ ^ 2 

It is not difficult to see tfiat inequalities (1.4) and (1.6) transmit to H as 
follows 

1 1 

6 ~6 



(6.8) |H(z,a)-H(z,6)Ki 



, for z G C and |^i|, |6| > R 



(6.9) sup|H(z,OI + supHflliLJi^!^ ^6M, for |^| > i? 
zen zi^z2 \z1-z2r 

The verification of (6.9) is a routine matter of the triangle inequality. 

The desired solutions : D — t- C of (6.1) will be obtained as restrictions 
to the unit disk of the solutions, still denoted by f = :C-^C, of the 
extended equation 

(6.10) f, = U{z,X + f,) 

The advantage of passing to the extended equation lies in the use of singular 
integrals in the entire plane. We represent the solution in the form of the 
Cauchy transform oi u = fz 



(6.11) 



f{z) = Cooiz) ^ - 



1 1 
+- 

z — T r. 



coir) 



We search for the density function a; in a Besov space ^a{C) C ^p{C) , 
p> ^ > 2. The density function will be supported in the double disk 2D . 



Recall the well known inequality [3, Theorem 4.3.13] 

(6-12) 

Thus, 

(6.13) \fiz)\ ^ CpMp, 

6.2. The Besov Space ,^a(C). Let < a < 1 be the Holder exponent 
in (1.5). Let us choose and fix for the rest of our paper the integrability 
exponent 

3 



\Crj{z)\ ^ Cp\z\^-p \\rj\\p , z e C, r? G =Sff (C) 

for z G D 



3.14) 



P 



> 3 



a 



The Besov space ^aiC-) consists of functions uj G ^^(C) which satisfy 



(6.15) 



UJ 



def ||a;(- + r) -a;(-)||p ^ 

a,p = \\oj\\p + sup r^. < 00 



We have a continuous imbedding ^a(C) C =5f°°(C) together with a uniform 
bound [1, Theorem 7.34(c)] or [31, p. 131]. 



(6.16) 
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6.3. Beurling-Ahlfors Transform. Complex derivative fz of the func- 
tion in (6.11) is expressed by a singular integral, known as the Beurling- 
Ahlfors transform 

(6.17) Mz) = Su;{z) [[ -^^^ where a; = /, 

We denote the norm of the operator S : ^p{C) ^p{C) by 5p > 1. 
Thus we have 

(6.18) ^ 5p , hence || 5w ^ cj ||o,^p 
This combined with (6.16) yields 

(6.19) ll/zlloo = ll'^'^lloo ^SpBp\\u\\a,p 

6.4. The structural parameter Aq. We are now ready to reveal the con- 
ditions on the complex parameter A; namely, 



(6.20) |A| ^ Ao = max < 



^yl6Sp ^ 4 condition (Ai) 

^y 81 SpL ^ -s/SSpL condition (A2) 

120 Sp Bp M condition (A3 ) 

3R condition (A4) 

32 CpL^ 32L condition (A5) 



These are not the optimal numerical values, but they are chosen for the 
clarity in subsequent computations. 

6.5. Solving the Extended Equation. The equation (6.10) is now equiv- 
alent to the integral equation 

(6.21) uj = H{z,X + Suj) 

for a density function lo = oj^{z) , which we shall find uniquely in the set 

(6.22) S = {w: ||a;||a,p ^ 60 • M} C ^g(C) C ^^(C) 

This is a closed subset of JSf^'(C) . We shall view *8 as a complete metric 
space with respect to -norm. First observe that for each a; G 03 we 

have 

(6.23) ||5a;||oo ^ GOSpEpM ^ ^Ao , 

by (6.19) and condition (A3). This combined with conditions (A3) and 
(A4) yields 

(6.24) \X + Soj{z)\ ^ ^|A| ^ ^\o>R 

In particular, the equation (6.21) is well defined on *8 . We now introduce 
a nonlinear operator T : *8 — > ^^(C) by the rule 

(6.25) Ta; = H(z,A + 5a;) 
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Clearly, suppTc<j C 2D . Note, by condition (6.9), that 

(6.26) WTojWp = \\U{z,X + Suj)\\p ^ 6M • ||x2d ||p = 6 (47r)"/^M ^ 18 • M 

Next we estimate the difference quotient for Tcj . To this effect we consider 
two cases 

Case 1 : ]r| ^ 1 . For all z G C, we can write 
\{Tu){z + t)-{Tu){z)\ 

(6.27) ^ |H(z + r, A + Suj{z + r)) - H(z + r, A + Su{z)) \ 
+ \H.{z + T,\ + Su{z)) - ^{z,\ + Suj{z))\ 

The first summand will be estimated by using (6.8) and (6.24), 

|H(z+r, \+Suj{z+t)) - H(z+r, \+Suj{z)) \ < ^ |<Sa;(^+r) - Suj{z) \ 

For the second summand in (6.27) we use Holder's estimate in (6.9) 

|H(z + r,A + 5a;(z)) - H(z, A + 5a;(z)) | ^ 6M|r|" 

We add both summands to obtain 

4L 

W 

The introduction of the factor X3d(^) legitimate beacuse the left hand 
side vanishes outside the disk 3D . We now compare the Jf^(C) -norms 
(with respect to z-variable) of both sides. 

II (Ta;)(- + r)-(Ta;)(-) ||p 

^ 1^ " ^^'^ + ■ 

Hence, by the condition (A2), 



|(Ta;)(^ + r) - {Tu){z) \ ^ ^ \Suj{z + r) - Suj{z) \ + 6M|rr • X^M 



(Tu;)(- + r)-(Ta;)(-)||p ^ 4L5p 



|r|" - |A| 

which together with (6.26) yields 



< — -rf • 60 M + 33 M ^ 36M 



(Tu;)(-+r)-(T.)(.)||, 

|r|° 

Case 2 : |r| ^ 1 . It suffices to use rough ^^-bound 



Irl 



11^^ ^3|M|,^54M ^60M 

Now we see that in both cases 1 and 2 we have a desired estimate || To; ||a^p ^ 
60 M, meaning that 

(6.28) T:<B^5S 

Next we show that T : — )■ *B is a contraction in the =5f^-norm. Let 
, a;2 G 03 . By (6.8) and (6.24) we see that 
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(6.29) 



I Twi — Tuj2 lip ^ L 

w 



A + Suji A + Slo2 
ll'5(^i-^2)|L 



Thus the contraction constant is at most ^ , because of the condition (A2). 
By virtue of Banach contraction principle the equation (6.21) has exactly 

one solution cj = G *B. We define f^{z) = Cuj^{z). 

It remains to examine the properties of /^. Using (6.16), (6.19) and (6.22) 
we estimate its derivatives as follows. 



and 



ll/g lloo ^ Bp \\co^\\a,p < 60BpM < ^Ao for ah z G C 



||/^1U = W^moo ^ SpBp \\ou%,p ^ mSpBpM ^ -Ao, z e 



Thus 
(6.30) 

which implies 
(6.31) 



^ ll/.^lloo + ll/i^lloo < Ao 



l/'(^)l < Ao, 



z G D 



This proves (6.3) in Proposition 6.1. Concerning continuity with respect to 
the parameter A, let |Ai|, IA2I ^ Ao- 

1 1 



00 



A2 



4L 



Ai + 5a;^i A2 + Suj^'^ 



ifP(2D) 



— ,l(Ai-A2) +5(a;^^-a;^^)||^,p^) 



4L 

16L 



Al - A2 ||^P(2D) + 



ALSp 



Al 



A2 



Al _ ^Aa 



Hence 



a; — a; 



A2 



^ IA1-A2I 



IA1A2 



Finally, we recall the estimate (6.12) for Cauchy transform, |C7y(2;)| ^ Cp \z 
for all z G C, which we apply to 77 = u^'^ — uj^'^ . This yields 



1- 



|/^^(.) - f>^\z)\ ^ ^|Ai-A2| ^Ao 



Al — A2 



Al A2 



, for all z G D 
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by the condition (A5), establishing the inequahty (6.4) in Proposition 6.1. 

□ 

We remark that is '^■'^'"-smooth. 

7. The difference of two solutions 

In this section we consider the difference of two solutions of (1.6), namely 
the given solution h and the good one ^'^(z) = Xz + f^{z)e constructed in 
the previous section. Precisely, let 

(7.1) g\z) = Xz + f\z) - h{z) 

We shall estimate the distortion of under the assumption that 

(7.2) ||/i,-(z)||^oo(D) =Ar<oo 

To this effect we must impose additional bound from below on the complex 
parameters A , which will now depend on the solution h . 

Proposition 7.1. Let h be as in Theorem 1.2 and as in Proposition 6.1. 
If 

(7.3) |A| ^4iV + 4Ao + i? 
then the difference function in (7.1) satisfies 

(7.4) I 55 (-2) I ^ 2 I ' almost everywhere in D 

Proof. Fix any point z G B at which and are defined. To simplify 
writing we omit the superscript A , so from now on f = f^ and g = g^ ■ 
Thus 

gz = fz - hz and g^ = X + f^ - 
We begin with a simple case 

Case 1: \hz\ ^ R- The computation goes as follows: 

2bK2(|M + |/i^-|) <2fjAo + iV 



^|A|-^Ao-ii^|A|-|M-|/i,| ^\g,\ 

Case 2: \hz\ > R, so we may apply the equation (1.6). First, we have the 
identity 

(A + fz) gz = {gz + hz) gz = gzgz+ K {fz - hz) 

Hence, by (6.6) 

(7.5) ^|A| • \gz\ ^ |A + fz\ ■ \gz\ sS \gz\ ■ \gz\ + \hz\ ■ \ fz - hz\ 

Here in the right hand side we estimate the term l^^^l by using (7.2) and (6.5) 

\g-z\ ^ \h-z\ + lA-K^+^Ao 
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For the second term in the right hand side of (7.5) we appeal to the equa- 
tion (1.6) and the Lipschitz condition (1.4), 

\h,\ ■ \h - h-,\ = \K\- \n{z, A + /,) - n{z, h,)\ 

For the last two estimates we have used (6.6) and condition (A2). Returning 
to (7.5), in view of |A + fz\ ^ ^|A| , we arrive at the desired estimate 

\X\-\g,\^{2N + 2Xo)-\g,\ 

That is, 

\gz\ < — ^ — - ■ \9z\ ^ 2 l^'^l ° 

Corollary 7.2. Under the assumptions of Proposition 7.1 the mapping 
is quasiregular, hence either constant or both open and discrete. 

7.1. Topological degree. Consider a continuous family of mappings 

(g^{z) = i g^{z) = z + j^ [f^{z) - h{z)] , for z G D and |A| ^ a 
jG°°(z) = z 

where 

(7.7) a = S\\h\\^oo(p) + 4||/i2||^oo(D) + 5Ao 
From (6.31) we have the uniform bound 

\G\z) - ^1 = ^ • |/'(^) - h{z)\ < ^ • (Ac + \\h\U < I 

By Corollary 7.2 the mapping G'^ : D — )■ C is quasiregular. Furthermore, 
it is nonconstant because for ^ < |2;| < 1 we have |G^(2;)| > — ^ > 
= |G^(0)|. Thus G^{z) is open and discrete. By Lemma 4.2 wc conclude 
that G^{zi) ^ G^{z2) , whenever zi and Z2 are distinct points in and 
|A| ^ (7 . This reads as follows 

Corollary 7.3. For all complex parameters A with |A| ^ a the mappings 
g^{z) = Xz + f^{z) — h{z) are injective in the disk |D ; that is, for z\ 7^ Z2 
in |D 

(7.8) h{z^) - h{z2) + X{zi - Z2) + fHzi) - f\z2) 

We shall infer from this, using topological degree arguments, the following 
inequality 



Lemma 7.4. For every circle = {A: |A| = p} with p ^ a there exists 
X ETp such that 

(7.9) \h{zi) - h{z2)\ ^ \X{zi-Z2) + f\zi)-f\z2)\ 
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Proof. This inequality certainly holds for large values of p . To simplify 
writing we denote a = h{zi) — h{z2) and assume, as we may, that a 7^ 0. 
We shall consider a family of mappings : T — > T , with parameter p ^ cr , 
given by 

= \Ftt^]-l\ ' ^^^'^ = A • (Zl - Z2) + f\zi) - f\z2) 

By virtue of the inequalities (7.8), each such mapping has well defined de- 
gree, denoted by deg<I>pi ^^Iso known as winding number. Letting the pa- 
rameter p vary wc obtain an integer- valued continuous function in p , thus 
constant. We identify this constant by letting p — > 00 . The mappings con- 
verge uniformly to : T ^ T , where $^(e*'') = j|^5f|| -e*^ . The degree 
of this limit map is equal to 1. Hence we conclude that 

deg $p = 1 , for all parameters p ^ cr 

We now fix p ^ C7 and move the point a 7^ to oc along the straight 
half-line {ta : t ^ 1 } , to observe that for some t ^ 1 the point ta lies 
in -F(Tp) . For if not, we would have well defined degree of the mappings 
: T ^ T , given by 



P ' ' |F(pe^^) - ta\ 
By virtue of continuity with respect to the parameter t we would have 
deg$*" = deg$^ = 1 , for alH ^ 1 

On the other hand letting t — 00 the mappings : T — > T converge 
uniformly to a constant map = , whose degree is zero, in contra- 
diction with the case t = 1 . Thus ta G F(Tp) , for some i ^ 1 , meaning 
that 

(7.10) ta = \-{zi- Z2) + f^{zi) - f^{z2) , for some A G Tp 

which yields the desired inequality (7.9). □ 

8. Proof of Theorem 3.1: the final step 
We invoke Lemma 7.4 with p = o", where a is given by (7.7). Thus 

IM^i)-M^2)Kt|a| ^^1^1-^21 + 1/"^'' (^i)-/"^'" (^2)1 ^ (a + Ao)|zi-Z2| 

where the latter inequality follows from Proposition 6.1. In terms of the 
gradient of h this reads as, 

(8.1) \Vh{z)\ ^a + Xo = 3\\h\y 00 (^D) + ^Wh^^oo^p^ + 6 Xo 

The final procedure, like for the Hopf-Laplace equation, consists of rescaling 
the variables. 

Let B(a, r) = {z : \z — a\ ^ r ^ 1} be a closed disk contained in il and 
h a solution to the equation hz = T-L{z,hz) in Consider the function 
h{z) = r~^h{rz + a) . It satisfies an equation hz = T-L{z, hz) in the unit 
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disk, where T-L(z,(,) = T-L{rz + a, ^) . We stress that < r ^ 1 , so aU the 
conditions imposed on Hiz,^), including (6.9), transmit to H(^;,^), with 
the same constants; easily verified. In other words the structural parameter 
Ao remains the same. Now inequality (8.1) applied to h yields 

3 

(8.2) l|V/i||_2»oo(iB) ^ ~ ll^ll.if°°(B) + 4 ||/lz||^oo(B) + 6Ao 

One may subtract a suitable constant from h, arriving at (3.2). □ 
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